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Abstract 

The charged-particle muhiphcity distribution is measured for all hadronic events 
as well as for light-quark and b-quark events produced in e"^e~ collisions at the Z 
pole. Moments of the charged-particle multiplicity distributions are calculated. The 
Hq moments of the multiplicity distributions are studied, and their quasi-oscillations 
as a function of the rank of the moment are investigated. 
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Introduction 



Since quarks and gluons are not observed directly, the understanding of the hadronization 
process whereby a quark-gluon system evolves to hadrons is of importance and provides a tool 
for studying the quark-gluon system itself. One of the most basic characteristics of the resulting 
hadronic system is the distribution of the number of hadrons produced. 

Assuming local parton-hadron duality (LPHD) P], characteristics of the charged-particle 
multiplicity distribution are directly related to the characteristics of the corresponding parton 
distributions. The parton distributions are calculable using perturbative quantum chromo- 
dynamics (pQCD). In particular, the dependence on the center-of-mass energy, -y/s, of the mean, 
(n), of the charged-particle multiplicity is an important test of pQCD. Since these calculations 
are only valid for light quarks, a separate measurement for light quarks is of interest. 

In this Letter, the charged-particle multiplicity distributions of hadronic decays of the Z 
boson are measured for b- and for light-quark (u, d, s and c) events as well as for all events. From 
these distributions moments are calculated, which characterize the shape of the distributions. 

The shape of the charged-particle multiplicity distribution is a fundamental tool in the 
study of particle production. Independent emission of single particles leads to a Poissonian 
multiplicity distribution. Deviations from this shape, therefore, reveal correlations 0. To study 
the shape, we use the normalized factorial moments. In terms of the multiplicity distribution, 
P(n), the normalized factorial moment of rank q is defined by 



It reflects correlations in the production of up to q particles. If the particle distribution is 
Poissonian, all Fq are equal to unity. If the particles are correlated, the distribution is broader 
and the Fq are greater than unity. If the particles are anti-correlated, the distribution is 
narrower and the Fq are less than unity. 

Normalized factorial cumulants, Kq, obtained from the normalized factorial moments by 



measure the genuine correlations between q particles, i.e., g-particle correlations which are not 
a consequence of correlations among fewer than q particles. 

Since \Kq\ and Fq both increase rapidly with g, it is useful to define the Hq moments. 



which have the same order of magnitude over a large range of q. 

The shape of the charged-particle multiplicity distribution analyzed in terms of the Hq was 
found to reveal quasi-oscillations [SlIllEllE], when plotted versus the rank g, in e"'"e~, as well 
as hadron-hadron, hadron-ion and ion-ion interactions. In e"'"e~ annihilation, this result was 
interpreted OHI in terms of pQCD, from which the Hq of the parton multiplicity distribution 
were calculated jSlEI- The expected behavior of Hq vs. q is quite sensitive to the approximation 
used, as is illustrated quahtatively in Figure [T] for the double logarithm approximation (DLA), 
the modified leading logarithm approximation (MLLA), the next-to- leading logarithm approx- 
imation (NLLA), and the next-to-next-to- leading logarithm approximation (NNLLA). In the 
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NNLLA a negative first minimum is expected near g = 5 and quasi-oscillations about zero are 
expected for larger values of q. 

According to the LPHD hypotliesis, liadronization does not distort the shape of the multi- 
phcity distribution. If this is vahd, the same shape may be expected for the charged-particle 
multiphcity distribution as for the parton multiphcity distribution. 

Experimental procedures 
Event selection 

This analysis is based on 1.5 million hadronic events collected by the L3 detector jH] at lep in 
the years 1994 and 1995 at the Z pole. 

Events are selected in a two-step procedure (TU]. First, at least 15 calorimetric clusters of 
at least 100 MeV are required in order to reduce background from the e"'"e~ t^t~ process. 
Hadronic events from the process e"'"e~ qq are then selected by requiring small energy 
imbalance both along and transverse to the beam direction. 

The second step is the selection of charged tracks measured in the central tracker and the 
silicon micro-vertex detector. A number of quality cuts are used to select well-measured tracks. 
Further, the thrust direction calculated from the charged tracks is required to lie within the full 
acceptance of the central tracker. No selection specifically rejects or selects tracks from long- 
lived neutral particles. The track selection efficiency, determined from Monte Carlo, is about 
75%. The resulting data sample corresponds to approximately one million selected hadronic 
events, and has a purity of about 99.8%. 

To correct for detector acceptances and inefficiencies, we make use of the jetset 7.4 ^l] 
parton shower Monte Carlo program, tuned using L3 data. Events are generated, passed through 
the L3 detector simulation program J,2j, and further subjected to time-dependent detector 
effects. Then they are reconstructed and the events and tracks are selected in the same way 
as the data. For systematic studies we also use events generated by ARIADNE 4.2 [12]. For 
comparisons with the data we use HERWIG 5.9 ^1] as well as jetset. 

To select b- and udsc-quark enhanced samples, we use the full three-dimensional information 
on tracks from the central tracker to calculate for each track the probability that it originated 
at the primary vertex [15J. We select b- and udsc-quark samples with purities of about 96% 
and 93% and efficiencies of about 38% and 96%, respectively. 

Unfolding 

The resulting multiplicity distributions are fully corrected for detector resolution using an 
iterative Bayesian unfolding method JH] • The detector and generator level Monte Carlo events 
are used to construct a matrix R{ndet, n) which represents the probability that ridet tracks would 
be detected if n charged particles were produced. A distribution, Po(^), is assumed for n. For 
this Pq) the distribution expected in the detector is -Po'^*(ndet) = ^„ -Rl'^det, ''^)-Po('^)- This is 
compared to the actual distribution of the raw data, and, making use of Bayes' theorem, an 
improved multiplicity distribution is calculated, which replaces Po(^) in the above expression. 
This process is repeated iteratively until satisfactory agreement between the expected and 
actual raw data distribution is found. In practice, this occurs after the second iteration if the 
JETSET multiplicity distribution is chosen as Po('^)- 
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In addition, corrections are made for efficiency and acceptance of the event selection, initial 
state radiation, and Kg and A decays. Furthermore, the distributions for the b- and udsc- 
enhanced samples are corrected for the purity of the flavor selection. 

The unfolding method gives |TH] an estimate of the covariance matrix of the unfolded dis- 
tribution. This matrix, combined with the uncertainties on the corrections mentioned above, is 
used to determine the uncertainties on the moments of the multiplicity distribution for the all- 
and udsc-flavor cases. When the statistics is too small, as in the b-fiavor case, the uncertainty 
on the estimate of the covariance matrix is large. In this case we use a Monte Carlo method. 
Many Monte Carlo variations of the raw data multiplicity distribution are made, choosing the 
number of events at each multiplicity from a Poisson distribution having as mean the observed 
number of events. These Monte Carlo distributions are then analysed in the same way as the 
data distribution. The uncertainty on a moment is determined from the spread in values of the 
moments of the Monte Carlo distributions. For the high-statistics cases, both methods agree. 

Systematic uncertainties 

The following sources of systematic uncertainty are investigated: 

Selection. The value of each cut used in the event selection is varied independently over a 
reasonable range around the default value and the resulting fully corrected distributions, 
together with their covariance matrices, determined, and from them the moments of the 
multiplicity distribution. For each multiplicity, as well as for each multiplicity moment, 
we assign a systematic uncertainty of half of the maximum difference between the new 
values. The same procedure is followed for the track selection and flavor tagging. For 
flavor tagging there is an additional contribution due to an uncertainty of 2.5% in the 
purity of the resulting sample, which accounts for the different response of the tagging 
algorithm to data and Monte Carlo. 

Monte Carlo uncertainties. The analysis is repeated using ARIADNE instead of jetset 
to determine the corrections and the unfolding matrix. The difference between the two 
results is taken as the systematic uncertainty. Further, the c- and b-quark fragmentation 
parameters, ec and eb, are varied. Also, the strangeness suppression parameter is varied 
by an amount consistent with the measured Kg production rate jT2j. In each case, half the 
difference between the results using the two parameter values is taken as the systematic 
uncertainty. 

Unfolding method. Three contributions are determined: First, ARIADNE is used to derive 
the initial distribution. Secondly, the analysis is repeated using a different number of 
iterations in the unfolding. Finally, the detector level multiplicity distribution of events 
generated by ARIADNE is unfolded using the response matrix, R{niet,'n), determined using 
JETSET events. In each case, the difference from the default value is taken as the systematic 
uncertainty. 

Background. The background of about 0.2% is mostly from two-photon processes. We take 
as a systematic uncertainty the effect of twice the amount of estimated background. 

The contributions from each of these sources are added in quadrature. The track selection 
contributes the dominant part of the total systematic uncertainty when all events are used, 
while the flavor-tagging purity uncertainty dominates that of the udsc sample. For the b-quark 
sample, these two contributions are about equal. 
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Source 


full 

lull 


udsc 


D 


Event selection 


u.uuo 


u.uuo 


n nn/i 

U.UU4 


Track selection 


n non 
u.uyu 


n nsn 
U.UoU 


U.iiO 


idgging cuts 




U.UiO 


n noi 

U.UZi 


Tagging purity 




U.ioO 


n 1 0R 

U.iZO 


MC modeling 


U.UC)/ 


U.Uoi 


n n/in 
U.U4U 


Unfolding 


U.Uo4 


U.Uo4 


U.U40 


RaclcPTonnd 


0.024 


0.024 


0.023 


7 conversion 


0.039 


0.039 


0.039 


Total 


0.11 


0.21 


0.19 



Table 1: Contribution of the various sources of systematic uncertainty to the measurement of 
the mean charged-particle multiplicity, (n). 

In addition, the accuracy of the simulation of the rate of photon conversion is considered. 
This is found to be about 15% smaller than in data jTHj and is assigned as a systematic 
uncertainty on (n). It is found to be negligible for the other moments. Breakdowns of the 
systematic uncertainties on {n) are shown in Tabled 



Results 

Charged-particle multiplicity distributions 

Charged-particle multiplicity distributions are measured both including and excluding Kg and 
A decay products.*^ Figure El shows the charged-particle multiplicity distribution including Kg 
and A decay products for the full, udsc- and b-quark samples. All distributions agree rather 
well with JETSET, but in all cases herwig gives a poor description of the data, as is seen in 
Figures EK andEb- 

From these distributions various moments of the charged-particle multiplicity distribution 
are calculated. The results are summarized in Table El The mean multiplicity including Kg and 
A decay products is consistent with our previous measurements [THl[Tn| and about 0.6 below 
the world average (21.07 ± 0.11) J7]. The difference in mean multiplicity between the cases of 
including or not the Kg and A decay products is consistent with our measurement of the K° and 
A production rates [2111 and with the world average ^12:. All the moments, with the exception 
of the dispersion, D, show significant flavor dependence. However, the flavor dependence of F2 
is quite small. F2 is also quite insensitive to the inclusion or not of Kg and A decay products. 
The difference between the mean charged-particle multiplicity of the b-quark sample and that 
of the udsc-quark sample is 2.58 ± 0.03 ± 0.08 when Kg and A decay products are included and 
2.43 ± 0.03 ± 0.08 otherwise. 



*''Note that S^, S^and fl^ have only one charged particle among their decay products apart from those 
produced in A decay, and E'^ and 'EP have none. Thus including or not the decay products of these baryons 
does not affect the charged multiplicity except through the A decay. 



4 



The Hq are calculated from the unfolded charged-particle multiplicity distributions. Since the 
Hq are sensitive to low statistics at very high multiplicities, we truncate the multiplicity distri- 
bution. The Hq thus obtained are biased estimators of the Hq of the untruncated distribution. 
This bias increases with stronger truncation, while the statistical uncertainty decreases, which 
allows a more significant comparison with models. It was suggested that even without 
this truncation, the Hq may be biased since a natural truncation occurs as a consequence of 
the finiteness of the sample. The truncation can induce oscillations or increase their size |21j . 
The truncation also introduces correlations between the Hq, although these are small for low 
q pT|[TUll^ . We choose the point of truncation such that multiplicities with relative error on 
P{n) greater than 50% are rejected. This corresponds, for all multiplicity distributions studied, 
to about 0.005% of events. For all three samples (full, udsc, and b) the truncation is at 53 if 
Kg and A decay products are included in the multiplicity and at 49 when they are not. The Hq 
presented here are calculated from distributions not including these decay products. However, 
the Hq are insensitive to their inclusion [TUj . 

The Hq of the truncated charged-particle multiplicity distribution from all, udsc- and b- 
quark events, shown in Figure El have a first negative minimum at g = 5 and quasi-oscillations 
for larger q. They are very similar for the three samples, with only slight differences for the 
b-quark sample. Similar behavior is seen for jetset (Figure Et)- Oscillations are also observed 
for HERWIG (Figure Eli), but they do not agree with those seen in the data. For both data and 
the Monte Carlo models, truncation at a lower value increases the depth of the first minimum 
and the amplitudes of the oscillations, while truncation at a larger value has the opposite effect. 

We note that our Hq, based on an order of magnitude greater statistics, agree with the Hq 
of SLD if we truncate at a value equal to the maximum multiplicity they observed |^. 

No truncation, other than that due to the finiteness of the sample, reduces the amplitudes 
of the oscillations to statistical insignificance, but the minimum at g = 5 remains, as is shown 
in Figure El Again, jetset agrees well with the data, while HERWiG does not. 

To investigate the effect of sample size on the Hq, 100 samples of pythia [SB] Monte Carlo 
events, were generated for sample sizes of 10^, 10^ and 10^ events, and their Hq determined. 
Their ±1 standard deviation bands are shown in Figure 01 In the insert of Figured the mean 
of the values is shown. For large q the values of the Hq depend on the sample size. However, 
for small q the values of the Hq are stable. In particular, H^ (the first minimum) changes little 
with the sample size, giving us confidence that the measured H^ is robust. Figure IH suggests 
that at least 10^ events, an order of magnitude beyond the statistics of the present experiment, 
would be needed to establish the maximum at g = 8. 

Conclusions 

The charged particle multiplicity distribution of hadronic Z decay and its moments are measured 
for light-quark and for b-quark, as well as for all flavor events. The Hq moments of truncated 
multiplicity distributions, which have smaller statistical uncertainties than those of the full 
distributions, are plotted versus the rank q. A negative minimum is observed at g = 5 followed 
by quasi-oscillations about zero, which is qualitatively similar to the behavior expected in 
NNLLA for the Hq moments of the full multiplicity distribution. Since Monte Carlo studies 
show that these oscillations are magnified, or even created, by truncation of the multiplicity 
distribution, the Hq are also measured for the untruncated multiplicity distribution. In this 



5 



case the minimum at g = 5, expected in both MLLA and NNLLA, is confirmed. But the 
oscillations at higher values of g, which are expected only in NNLLA, cannot be confirmed. 
Previous observations of these oscillations are most likely a consequence of truncation resulting 
from limited statistics. 
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All events 


without Kg and A decay 


with Kg and A decay 


(n) 
in') 

(n^) . 10-2 

{n") ■ 10-3 

D=^{{n-{n)y) 

S = {{n- {n)f)/D^ 

K^{{n-{n)Y)/D^-?, 

{n)/D 

F2 = {n{n - l))/(n)2 


18.63 ± 0.01 ± 0.11 
381.7 ± 0.3 ± 4.4 

85.2 ± 0.1 ± 1.5 
205.9 ± 0.4 ± 5.1 
5.888 ± 0.005 ± 0.051 
0.596 ± 0.004 d= 0.010 

0.51 ± 0.01 ± 0.04 
3.164 ± 0.002 ± 0.016 
1.0461±0.0002± 0.0040 


20.46 ± 0.01 ± 0.11 
457.7 ± 0.3 ± 4.9 
111.1 ± 0.1 ± 1.8 
290.6 ± 0.5 ± 6.5 
6.244 ± 0.005 ±0.051 
0.600 ± 0.004 ±0.010 
0.49 ± 0.01 ± 0.03 
3.277 ± 0.002 ±0.016 
1.0441±0.0001±0.0034 


udsc-quark events 






(n) 
in') 

(n^) ■ 10-2 

(n^) ■ 10-3 

D= v/((n-(n))2) 

S= {{n-{n)f)/D^ 

K= {{n- {n)y)/D^-3 

{n)/D 

F, = {n{n-l))/{ny' 


18.07 ± 0.01 ± 0.21 
340.0 ± 0.3 ± 8.4 

78.3 ± 0.1 ± 2.7 
184.4 ± 0.4 ± 8.6 
5.769 ± 0.007 ± 0.071 
0.613 ± 0.005 ± 0.014 

0.54 ± 0.02 ± 0.06 
3.133 ± 0.003 ± 0.020 
1.0161±0.0002± 0.0015 


19.88 ± 0.01 ± 0.21 
432.4 ± 0.4 ± 9.2 
102.2 ± 0.1 ± 3.3 
260.7 ± 0.5 ± 11.1 

6.111 ± 0.007 ± 0.071 
0.617 ± 0.005 ±0.012 
0.53 ± 0.02 ± 0.05 
3.252 ± 0.003 ± 0.020 

1.0111±0.0002±0.0038 


b-quark events 






(n) 
(n2) 

(n3) . 10-2 
(n^) • 10-3 
D^^{{n-{n)y) 
S= {{n-{n)f)/D^ 
K= {{n- {n)y)/D^-3 
{n)/D 

F2 = {n{n - l))/(n)2 


20.51 ± 0.02 ± 0.19 
453.9 ± 1.1 ± 1.8 
107.9 ± 0.4 ± 0.7 
273.9 ± 1.5 ± 1.9 

5.78 ± 0.01 ± 0.07 
0.574 ± 0.017 ± 0.008 

0.43 ± 0.04 ± 0.04 
3.551 ± 0.006 ± 0.055 
1.0305±0.0003± 0.0027 


22.45 ± 0.03 ± 0.19 

542.0 ± 1.2 ± 3.0 

140.1 ± 0.5 ± 1.1 
385.8 ± 1.9 ± 1.7 

6.16 ± 0.01 ± 0.07 
0.573 ±0.017 ±0.007 

0.42 ± 0.04 ± 0.03 
3.645 ± 0.005 ± 0.049 
1.0307±0.0002±0.0023 



Table 2: Moments of the charged-particle multiplicity distribution for all, udsc-, and b-quark 
events. The first uncertainty is statistical, the second systematic. 
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Figure 1: Qualitative behavior of Hq as a function of q for various approximations 
of perturbative QCD [3,8]. 
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Figure 2: Charged-particle multiplicity distribution for all, udsc-, and b-quark 
events compared to the expectations of (a,c,e,g) jetset and (b,d,f,h) HERWIG. The 
error bars include both statistical and systematic uncertainties. 
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Figure 3: The Hg of the truncated (a,b) and non-truncated (c,d) charged-particle 
multiphcity distribution compared to the expectations of (a,c) jetset and (b,d) 
HERWIG. The error bars include both statistical and systematic uncertainties. 
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Figure 4: 1-standard deviation bands of expectated Hq of the non-truncated 
charged-particle multiplicity distribution from pythia for sample sizes of 10^, 10^ 
and 10^. The insert shows the mean Hq of 100 samples of 10^, 10^ and 10'' pythia 
events. 
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